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ABSTRACT 


Bourdon  tubes  are  the  most  commonly  used  elastic  elements  in  mechanical 
pressure  gauges.  Basic  theory  describing  the  principles  of  Bourdon  behavior  is  readily 
available  but  very  few  analytical  studies  have  been  published  that  model  Bourdon 
element  behavior.  The  purpose  of  this  research  is  to  develop  an  analytical  model  to 
determine  tip  displacement  in  tubular  elastic  elements  using  basic  principles  of  solid 
mechanics.  To  determine  the  validity  of  the  results  obtained  from  the  analytical  model,  a 
study  of  results  from  finite  element  solutions  and  experimental  data  is  also  presented. 
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I.  INTRODUCTION 


A.  MOTIVATIONS  FOR  TfflS  RESEARCH 

Bourdon  tubes,  the  elastic  elements  used  in  mechanical  pressure  gauges,  are 
quite  popular,  even  today,  due  to  their  simplicity.  When  a  Bourdon  tube  is  fabricated 
from  conventional  alloys  and  operated  within  its  elastic  limit,  the  tip  of  the  tube  moves  a 
small  distance.  However,  if  the  Bourdon  element  were  fabricated  from  a  superelastic 
alloy  capable  of  withstanding  strains  on  the  order  of  8%,  strains  much  higher  than  what 
ordinary  alloys  can  withstand,  the  Bourdon  tube  could  be  used  in  unconventional 
applications  as  a  robotic  finger  or  as  an  actuator.  With  these  imique  applications  in 
mind,  the  study  of  the  elastic  deformation  of  a  Bourdon  tube  was  initiated. 

The  efforts  of  a  literature  survey,  the  results  of  which  are  discussed  later  in  this 
chapter,  yielded  little  information  in  terms  of  analytical  solutions  to  Bourdon  tube 
behavior.  In  industries  that  manufacture  Bourdon  tubes,  tip  travel  is  obtained  from 
experimental  data  in  the  absence  of  a  mathematical  model.  This  further  prompted  the 
development  of  an  analytical  solution  for  tip  travel  in  a  Bourdon  tube. 

With  today’s  micro-chip  technology,  many  people  think  that  the  trend  for  pressure 
gauges  is  moving  in  the  direction  of  the  electronic  transducer.  Surprisingly,  precision 
pressure  gauges  demonstrate  consistent  and  accurate  behavior  using  Bourdon  tube 
technology.  As  a  result,  the  demand  for  Bourdon  gauges  is  quite  significant.  The  most 
obvious  advantages  characteristic  of  Bourdon  gauges  are  their  (1)  broad  pressure  range, 
(2)  low  cost  compared  to  electronic  gauges,  (3)  ease  of  manufacturing,  (4)  ease  of 
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operation,  (5)  ease  of  troubleshooting,  (6)  good  accuracy  (up  to  0.1  %  full  scale 
deflection),  and  (7)  lack  of  an  external  power  source  which  precludes  susceptibility  to 
line  noise  and  power  outages.  [Ref  1,2] 

Bourdon  tubes  are  extremely  popular  as  elastic  elements  in  pressure  gauges. 
Despite  their  popularity  little  analytical  work  has  been  performed  or  published  modeling 
their  behavior.  The  work  presented  in  this  thesis  provides  such  an  analysis. 


B.  BOURDON  TUBE  BEHAVIOR 

1.  Developmental  History 

In  1849,  a  French  engineer  named  Eugene  Bourdon  invented  the  C-shaped 
element  commonly  known  today  as  the  Bourdon  tube  or  C  spring.  While  witnessing  the 
construction  of  a  steam  engine,  Eugene  Bourdon  noticed  that  the  helically  wound  coil 
used  to  condense  the  steam  was  flattened  dming  formation.  To  correct  the  deformation, 
the  tube  was  plugged  at  one  end  and  pressurized  at  the  other.  As  a  result,  the  coil 
partially  unwound  as  the  flattened  tube  regained  its  circular  cross-section.  Intrigued  by 
what  he  saw,  Eugene  Bourdon  conducted  his  own  experiments  and  eventually  invented  a 
pressure  gauge  using  the  tip  deflection  of  a  curved  tube  with  an  elliptical  cross  section  as 
a  pressure  measuring  device.  [Ref  1]  Figure  1  offers  an  exploded  view  of  a  typical  C- 
type  pressure  gauge. 

2.  Deflection  Theory 

The  Bourdon  tube,  henceforth  referred  to  as  the  Bourdon,  is  the  actuating  element 
in  many  pressure  measuring  devices.  Typically,  the  Bourdon  elastic  element  has  an  oval 
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or  elliptical  cross-section  and  is  coiled  in  the  shape  of  a  C,  a  spiral,  or  a  helix  as  shown  in 
Figure  2.  One  end  of  the  tube  is  sealed  and  the  other  is  connected  to  a  capillary  tube 
through  which  pressure  is  applied.  Any  pressure  within  the  tube  in  excess  of  the  external 
or  atmospheric  pressure  causes  the  flattened  element  to  elastically  change  its  shape  to  a 
more  circular  cross-section.  The  resulting  change  in  the  cross-section  causes  the  coiled 
tube  to  straighten.  With  one  end  of  the  tube  fixed,  the  straightening  effect  causes  the  free 
end  to  move  a  small  distance.  This  small,  imenhanced  movement  is  called  tip  travel  as 
shown  in  Figure  3.  [Ref  1,3,4,5,6] 

In  a  C-type  Bourdon,  the  elastic  deformation  of  the  cross-section  creates 
compressive  stresses  along  the  inner  wall  and  tensile  stresses  along  the  outer  wall  as 
defined  in  Figure  3.  When  added  to  the  tensile  stresses  existing  in  both  walls  due  to  the 
applied  internal  pressure,  a  stress  imbalance  is  created  such  that  the  tensile  stress  in  the 
outer  wall  is  greater  than  the  tensile  stress  in  the  inner  wall.  In  the  process  of 
straightening,  the  tube  experiences  tensile  stresses  in  the  inner  wall  and  compressive 
stresses  in  the  outer  wall,  just  the  opposite  caused  by  the  change  in  cross-section.  The 
effect  of  straightening  relieves  the  unbalanced  stress  condition  created  by  the  change  in 
cross-section.  The  free  end  of  the  Bourdon  moves  with  the  application  of  pressure  until 
opposing  stresses  within  the  tube  reach  equilibrium.  [Ref  1,4,5] 

3.  Design  Criteria 

As  stated  earlier,  the  first  Bourdon  pressure  gauge  was  patented  more  than  100 
years  ago.  But  despite  its  age  and  the  fact  that  “technology”  is  now  much  more 
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Figure  3.  Description  of  tip  travel  in  a  Bourdon  element  (from  Ref  1) 
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advanced,  very  little  about  the  design  of  a  Bourdon  gauge  has  changed.  The  factors 
considered  most  important  in  the  design  of  tubular  elastic  elements  are  (1)  the  element 
material,  (2)  the  coil  length,  (3)  the  wall  thickness  of  the  tube,  (4)  the  cross-sectional 
shape,  and  (5)  the  desired  tip  travel. 

Bourdons  may  be  made  of  ferrous  or  nonferrous  metals  depending  on  the 
application  of  the  gauge.  Admiralty  brass,  phosphor  bronze,  and  beryllium  copper,  which 
are  used  to  make  grade  B  accuracy  gauges,  higher  quality  grade  A  and  AA  gauges,  and 
high  pressure  gauges,  respectively,  are  examples  of  nonferrous  metals  applications.  Ni- 
Span  C  is  a  popular  but  expensive  ferrous  material  used  for  its  exceptional  temperature 
behavior.  [Ref  1] 

As  stated  earlier.  Bourdon  elements  may  be  coiled  in  the  shape  of  a  C,  a  spiral,  or 
a  helix.  The  shape  and  coil  length  of  the  elastic  element  depends  on  the  pressure  range 
of  operation.  For  pressures  up  to  60  bar  (6  MPa,  870.24  psi),  C-type  elements  are 
typically  used  and  for  pressures  above  60  bar,  spiral  or  helical  elements  are  used. 

[Ref  5] 

The  geometry  of  the  cross-section,  one  factor  affecting  tip  displacement,  is 
chosen  depending  on  the  pressure  range  of  operation.  The  distinction  is  made  between 
low  and  high  pressure  measuring  elements  for  the  proper  choice  of  the  wall  thickness. 
Low  pressure  gauges  are  loosely  defined  as  those  operating  from  12  -  500  psi  and 
typically  have  flatter  cross-sections  and  thirmer  walls  than  high  pressure  gauges  which 
have  rounder  cross-sections  and  thicker  walls  and  may  read  pressures  as  high  as  100,000 
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psi.  [Ref.  6]  Low  pressure  Bourdons  experience  greater  tip  travel  for  a  given  pressure 
whereas  high  pressure  Bourdons  exert  higher  tip  forces.  [Ref  1] 

Tip  travel  in  a  Bourdon  element  varies  widely  with  the  range  of  the  gauge.  Low 
pressure  gauges,  with  thinner  walls  and  smaller  minor  axes,  are  less  effective  in 
overcoming  the  frictional  forces  within  the  amplifying  movement.  Higher  tip  travel  and 
smaller  amplification  means  smoother  pointer  motion  but  also  produces  excessive 
stresses  in  high  pressure  elements.  As  a  result,  high  pressure  gauges  have  thicker  walls, 
larger  minor  axes,  and  high  amplification  movements  to  keep  tip  travel  small.  Gauge 
manufacturers  design  their  Bourdon  elements  for  specific  magnitudes  of  tip  deflection. 
For  ease  of  manufacturing  and  maintenance,  variations  in  design  are  kept  to  a  minimum. 
As  a  result,  the  major  distinction  between  gauges  is  found  within  the  movement 
mechanisms. 

When  designing  their  gauges,  manufacturers  either  determine  the  appropriate 
length  of  the  sector  tail  and  the  gear  ratio  necessary  to  produce  the  pointer  rotation  for 
the  desired  tip  travel  or  determine  the  necessary  tip  travel  to  operate  a  specific 
movement.  [Ref  1,2,  3,  4,  5]  Bourdon  tip  displacement  is  generally  linear  with  pressure 
such  that  25%,  50%,  and  75%  of  the  total  applied  pressure  is  equal  to  25%,  50%,  and 
75%  of  the  total  tip  displacement.  Non-linearities  in  tip  displacement  are  compensated 
by  varying  bourdon  geometry  or  by  adjusting  the  movement. 

4.  Temperature  Effects 

One  of  the  most  difficult  problems  that  manufacturers  of  precision  pressure 
gauges  must  overcome  is  the  effects  of  temperature  on  the  elastic  element.  Tip 
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displacement  is  very  susceptible  to  changes  in  ambient  temperature.  The  temperature 
characteristics  of  a  Bourdon  depend  on  geometry  and  material.  Fluctuations  in 
temperature  affect  the  elastic  modulus,  or  stiffness,  of  the  bourdon.  Higher  temperatures 
correspond  to  lower  stiffnesses  which  mean  greater  tip  travel.  The  opposite  is  also  true. 

Thermal  length  changes  affect  the  entire  measuring  system.  The  easiest  and 
cheapest  way  to  compensate  for  temperature  effects  is  to  manufacture  elastic  elements 
from  materials  having  low  coefficients  of  thermal  expansion.  The  material  most 
commonly  used  for  this  purpose  is  Ni-Span  C.  Although  less  affected  by  temperature, 
Ni-Span  C  does  not  completely  eliminate  the  effects  of  temperature.  Manufacturers  also 
compensate  for  temperature  effects  by  making  all  of  the  components  of  the  measuring 
device  from  the  same  material.  This  approach  eliminates  any  tip  deflection  relative  to 
the  movement  caused  by  temperature  changes.  [Ref  2,4] 

Some  gauge  manufacturers  maintain  gauge  accuracy  across  fluctuating 
temperatures  with  the  addition  of  a  temperature  compensator  within  the  movement 
mechanism.  The  device  reacts  to  the  same  temperatme  that  affects  the  bourdon.  As  the 
temperature  changes,  the  device  adjusts  a  setting  within  the  movement  to  offset  the 
temperature  effect.  Each  compensator  is  custom  manufactured  to  match  and  offset  the 
temperature  characteristics  of  the  elastic  element  with  which  it  is  used.  [Ref  2] 

5.  Technical  Publications 

Only  one  publication  was  found  that  presented  an  analytical  solution  for  bourdon 
behavior.  In  1946,  Alfred  Wolf  wrote  An  Elementary  Theory  of  the  Bourdon  Gage  which 
presented  theory  based  on  two  elements  of  strain;  strain  due  to  “the  bending  of  the  walls 


9 


in  a  transverse  section  through  the  tubing”  and  to  “longitudinal  extension  and  contraction 
parallel  to  the  axis  of  the  tubing”.  In  this  publication.  Wolf  derived  formulas  for  what  he 
called  the  sensitivity  {SR/R)  and  the  torque  (7)  of  the  bourdon  where  SR  was  defined  as 
the  change  in  radius  of  curvature  and  torque  was  defined  as  the  external  bending  moment 
required  to  restore  tip  deflection  to  zero  for  any  given  pressure.  Based  upon  his 
theoretical  results.  Wolf  concluded  that  the  theory  of  pure  bending  in  thin  shells  did  not 
adequately  model  the  behavior  of  the  bourdon  for  the  case  much  greater  than 

unity  and  that  gauge  sensitivity  was  inversely  proportional  to  the  first  power  of  t  for  thin- 
walled  tubes.  Wolf  s  tabulated  results  for  “observed”  and  “calculated”  sensitivities  and 
torques  for  a  steel  elastic  element  were  quite  comparable  to  on  another. 

6.  Manufacturers’  Literature 

The  logical  progression  in  the  search  for  literature  was  contacting  companies  that 
manufacture  the  elastic  elements  addressed  and  analyzed  in  this  paper.  A  representative 
from  the  Stewart  Warner  Instrument  Corporation  commented  that  their  gauges  cover  a 
very  broad  range  of  pressures  and  are  “old  standards”  based  on  tried  and  true  designs. 
On  the  rare  occasion  that  a  customer  requests  a  pressure  gauge  not  previously 
manufactured,  the  company  has  a  “formula”  for  estimating  the  new  design  based  on 
established  gauge  designs.  The  company  representative  also  said  he  had  no  knowledge 
of  an  analytical  model  to  determine  tip  travel. 

Excerpts  from  a  pressure  measurement  handbook  produced  by  the  WIKA 
Instrument  Corporation  provided  a  brief  history,  the  principles  of  operation,  some  design 
considerations,  the  effects  of  temperature,  and  load  constraints  of  bourdon  gauges.  The 
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excerpts  also  presented  equations  for  tip  travel,  relative  change  in  curvature,  and 
maximum  bending  stresses  of  tubular  elements  but,  for  proprietary  reasons,  all  were 
expressed  as  functions  of  undefined  variables. 

The  Heise  Division  of  Dresser  Industries  also  provided  some  literature  on  dial 
gauges.  Their  biggest  contribution  to  this  literature  survey  was  information  concerning 
temperature  compensation,  methods  of  which  were  previously  discussed.  Like  the 
Stewart  Warner  Corporation,  Heise’s  gauges  are  well  established  and  rooted  in  50-year 
old  designs.  Heise  is,  however,  currently  developing  its  own  analytical  model  for 
Bourdon  elements. 

The  efforts  of  a  literature  survey  yielded  relatively  little  information  concerning 
analytical  studies  of  Bourdon  pressure  gauges.  Very  few  technical  papers  have  been 
published  which  analytically  describe  Bourdon  theory.  As  a  result,  gauge  manufacturers 
typically  rely  on  experimentation  and  established  designs  when  creating  new  gauges. 


C.  RESEARCH  OBJECTIVES 

The  objectives  of  this  research  were  to  develop  an  analytical  solution  for  tip 
travel,  specifically  in  C-type  Bourdons,  and  to  compare  results  from  the  analytical  model 
to  results  from  finite  element  models  and  data  collected  from  laboratory  experiments. 
Because  pressure  gauge  manufacturers  most  commonly  rely  on  experimentation  for  the 
design  of  new  gauges,  a  realization  of  the  stated  objectives  could  prove  to  be  quite  useful 
to  them  or  to  anyone  interested  in  the  applications  of  Bourdon  theory. 
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n.  ANALYTICAL  ANALYSIS  OF  A  BOURDON  TUBE 


A.  NOMENCLATURE  AND  ASSUMPTIONS 
1.  Nomenclature 

For  clarity,  symbols  related  to  Bourdon  geometry  are  annotated  in  Figure  4  and 
defined  as  follows: 

a  -  semi-major  axis  of  the  elliptical  cross-section 
b  -  semi-minor  axis  of  the  elliptical  cross-section 
t  -  wall  thickness  of  the  Bourdon  element 
R  -  radius  of  curvature  of  the  Bourdon  element 
//  -  angle  describing  the  length  of  the  Bourdon  tube 
/  -  circumference  of  the  Bourdon  tube 
p  -  radius  of  curvature  of  the  cross-section  of  the  Bourdon  tube 
^  -  angle  describing  the  position  of  a  point  on  the  circumference  of  the 
Bourdon  element 

k  -  complete  modulus  of  elliptic  functions 
K  -  complete  elliptic  integral  of  the  first  kind 
E  -  complete  elliptic  integral  of  the  second  kind 
E‘  -  Young’s  modulus 

E,  -  elastic  modulus  from  the  theory  of  thin  shells 
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AZ 


Figure  4.  Bourdon  length  segment  and  cross-sectional  view 
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2.  Assumptions 

In  the  process  of  developing  this  analytical  model  for  Bourdon  behavior,  the 
following  general  assumptions  were  made: 

a)  the  circumference  /  remains  constant  during  cross-sectional  deformation, 

b)  the  length  of  the  tube,  defined  as  /jR,  remains  constant  as  the  elastic  element 
straightens, 

c)  the  wall  thickness  t  remains  constant  as  the  cross-section  changes, 

d)  the  change  in  the  angle  (j)  (^  -  ^)  is  so  small  for  small  deformations  that  ^ 
essentially  remains  constant  during  deformation,  and 

e)  the  material  of  the  elastic  element  is  homogeneous  and  isotropic. 

B.  MATHEMATICAL  BACKGROUND 

As  a  preface  to  this  section,  detailed  derivations  of  all  equations  are  presented  in 
Appendices  A-D.  For  the  sake  of  clarity  and  brevity,  only  simplified  expressions  are 
presented  in  the  body  of  this  paper. 

Some  background  concerning  Jacobian  elliptic  functions  and  integrals  is 
necessary  for  following  the  mathematical  progression  in  this  thesis.  The  following 
definitions 
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and 


*71 

2 


E= 


•Jo 


1  -  k^  sin-<j) 


2.2 


are  complete  elliptic  integrals  of  the  first  and  second  kind,  respectively,  [Ref  9,  10]  and 
appear  frequently  throughout  the  problem  formulation.  The  variable  k ,  defined  by  the 
expression 


2 

a 

on  the  interval  0<^<1,  is  known  as  the  modulus  of  the  Jacobian  elliptic  functions. 
Identities  containing  the  elliptic  functions  snfuj,  cn(u),  and  dn(u)  are  presented  in 
equations  2.4. a  -  2.4.f  which  reduce  to  the  trigonometric  functions  sinfu),  cos(u),  and  1 
when  k  is  zero.  Elliptic  function  identities  relevant  to  this  problem  formulation  are: 


sn^(u)  +  cn^(u)  =  l  2.4.a 

dn^(u)  +  k^sn^(u)=l  2.4.b 

sn(0)=0  2.4.C 

cii(0)=l  2.4. d 
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<in(0)=l 


2.4.e 


cn(K)=0  2.4.f 

which  are  also  defined  in  Refs.  9  and  10.  Solutions  for  the  following  specific  integrals 
are  derived  here  for  later  use  in  the  computation  of  strain  energy  in  elastic  Bourdon 
elements  due  to  cross-sectional  deformation; 

•K 

dn*'^  (u)  du* - - 

-0  3(1 -k') 


2E(2-k^) 

(1-k^) 


■K 

-0 


dn" 


*  (u)  du= 


4  (2 -k^) 
15(1 -kV 


2E(2-k^) 
(1  -  k") 


K 


3E 

5(1 -kV 


2.6 


■K 

dn'®  (u)  du= 

”0 


24-  (2  -  k^)^ 
105  (1  -  k^)^ 


5 

21- (1  -  kV 


2E(2-k^)  ^  18E(2-k^) 

(1-k^)  35(1 -kV 


2.7 
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C.  PROBLEM  FORMULATION 


1.  Geometric  Relationships 

Assuming  that  the  volume  of  the  Bourdon  changes  as  deformation  occurs,  the 
change  in  the  semi-major  axis  da  is  expressed  as  a  function  of  the  change  in  the  semi¬ 
minor  axis  db  by  the  equation 


ab(E-K) 

Ea^-Kb^ 


•db 


2.8 


The  derivation  of  equation  2.8  is  presented  in  Appendix  A.  For  a  unit  length,  the  volume 
of  the  Bourdon  is  defined  by  the  equation 


V*=7iab 


2.9 


The  change  in  volume  then  becomes 

dVssu(adb  +  bda) 


2.10 


or 


dV=7ia 


a^  +  b^)-E  -  2Kb^ 
Ea^  -  Kb^ 


db 


2.11 


by  making  the  substitution  for  da  in  equation  2. 10. 

Bourdon  tip  travel  is  a  function  of  the  sectional  deformation  that  occurs  as  a 
result  of  pressurization.  For  an  applied  internal  pressure,  a  curved  circular  tube  will  not 
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straighten.  But  a  Bourdon  straightens  as  its  elliptical  cross-section  becomes  more 
circular.  Because  of  this  behavior,  it  is  important  to  understand  how  the  radius  of 
curvature  of  the  cross-section,  p,  changes.  As  can  be  seen  in  Figure  4,  the  radius  of 
curvature  for  a  circular  section  is  constant  but  changes  as  a  function  of  the  angle  ^  for  an 
elliptical  section.  A  simplifying  assumption  here  is  that,  as  a  point  is  mapped  from  the 
elliptical  shape  to  the  more  circular  shape,  the  defining  angle  ^  is  approximately  the 
same  before  and  after  (^)  deformation  for  small  tip  travel  and  minimal  sectional 
deformation. 

Whereas  p  defines  the  radius  of  curvature,  1/p  defines  the  curvature  of  the  cross- 
section.  The  simplified  expression  for  curvature  may  be  written  as 


b 

a^  (i  -  k^  sirf  <t))*^ 


2.12 


The  change  in  curvature,  determined  by  taking  the  derivative  of  1/p  only  with  respect  to 
dimensions  a  and  h,  is  defined  by  the  expression 


W  a^  (l  -  k^  sin^  ';^)  (Ea^-Kb^) 

2.13 

The  derivative  of  1/p  was  not  taken  with  respect  to  the  angle  ^  because,  as  stated  earlier 
in  the  assumptions,  <ff  is  assumed  to  remain  constant  throughout  pressurization  and 
subsequent  deformation. 


C-a^-b^)E 


-K 


(1  -  K  sm  ■(/>)  +  S  E 


db 
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2.  Incremental  Strain  Energy  Due  to  Cross-Sectional  Deformation 

Having  derived  the  fundamental  geometric  relationships  necessary  to  develop  the 
analytical  model,  energy  methods  were  applied  to  obtain  a  solution  to  the  problem.  In 
general,  the  elastic  strain  energy  f/in  a  beam  due  to  bending  is  written  as 


rL 


u= 


-‘0 


M" 

2EI 

1 


dx 


2.14 


as  defined  in  Ref  1 1 .  For  an  initially  curved  beam,  the  bending  moment  A/ is  defined  as 


M=EIdW 

^  W 


2.15 


Figure  4  shows  that,  for  an  incremental  slice  of  the  cross-section,  the  wall  of  the  tube  can 
be  modeled  as  a  curved  beam  where  d(l/p)  is  the  change  in  curvature  of  the  cross- 
section.  From  the  theory  of  thin  shells,  the  effective  elastic  modulus  Ei  is  defined  by  the 
expression 


2.16 


where  E‘  is  the  elastic  modulus  of  a  straight  beam.  The  area  moment  of  inertia  1  is 
defined  for  a  unit  length  slice  of  the  tube  by  the  expression 


I=-  2.17 

12 
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Rewriting  the  strain  energy  equation  as  symmetric  in  four  quadrants,  the 
expression  for  the  total  incremental  strain  energy  due  to  cross-sectional  deformation  is 


dU  *4- 
xs 


2-EI 

1 


ds 


2.18 


Substituting  equations  2.15,  2.16,  and  A.  14  into  equation  2.18,  the  incremental  strain 
energy  becomes 


*71 

2 

Jo 


2.19 


and  eventually,  after  considerable  simplification. 


dU  * 


xs 


2aE^ICjdb^ 

•K 

C^dn  ‘^-(u)  +  C^dn  ^(u)  +  C^dn  *(u) 

(  2\ 

(1  -  V  )  J 

0 

du 


2.20 


where  the  constants  Cy  through  C4  are  defined  as  follows: 


a*(Ea^  -  Kb^)^ 


2.20.a 
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2 


2.20.b 


■E-K 


C  s6  E 

3 


\ 


•E-K 


/ 


2.20.C 


C  =9E^ 


2.20.d 


The  variables  AT  and  E  represent  the  complete  elliptic  integrals  of  the  first  and  second 
kind,  respectively,  and  are  defined  in  equations  2. 1  and  2.2.  The  integrals  of  (u), 
dn^(u),  and  dn'^(u)  are  defined  by  equations  2.5, 2.6,  and  2.7,  respectively. 

3  Incremental  Strain  Energy  Due  to  Bending 

To  determine  Bourdon  incremental  strain  energy  due  to  bending,  or  straightening 
in  this  case,  an  expression  was  developed  for  the  total  bending  moment  due  to  an  applied 
internal  pressure.  Referring  to  Figure  5,  the  net  moment  M  is  found  by  summing  the 
moments  acting  on  the  end  cap  and  walls  of  the  tube  about  point  O.  The  details  of  this 
derivation  are  shown  in  Appendix  C  and  the  resulting  expression  is 

— ^ 

M=-P  a-R  D  (l  -  cos^)  k  2.21 

where  P  is  the  applied  pressure,  R  is  the  radius  of  curvature  measured  to  the  centerline  of 
the  tube,  and  £>  is  a  constant  for  specified  dimensions  a  and  b  defined  as 
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end  cap;  Bourdon  tip 


Figure  5.  Bourdon  length  segment  and  cross-sectional  view  depicting  the  elemental  area 
for  the  derivation  of  the  net  moment  acting  on  the  elastic  element 
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4 

K  -  E  ,  _ 

b' 

— — +  (2E-K) 

+  X  — 

3 

k' 

a 

Again,  K  and  E  represent  complete  elliptic  integrals  of  the  first  and  second  kind, 
respectively,  as  defined  in  equations  2.1  and  2.2.  Substituting  equations  2.21  into 
equation  2.14,  the  total  strain  energy  due  to  bending  of  a  length  of  tube  that  subtends  the 
angle  //  with  the  center  of  curvature  becomes 


(-P  a  R  D  (i  -  cos(p)))^ 
2E‘I 


Rdp 


-'0 

where  the  moment  of  inertia  /  is  now  defined  as 


I*~  ((a  +  t)  (b  + 1)^  -  a  b^) 


2.23 


2.24 


The  incremental  strain  energy  due  to  bending  was  foimd  by  differentiating  C4  with 
respect  to  a,  b,  P,  and  R  as  follows: 

dU.  =1  u.  +  i-  U,  +  ^  U.  +  —  U.  2.25 

b  b  b  b  b 

da  db  dP  dR 

To  simplify  the  results  of  this  lengthy  equation,  the  incremental  strain  energy  due  to 
bending  is  rewritten  as 
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dU^=A(B  +Z) 


2.26 


where  A,  B,  and  Z  are  defined  as  follows: 


A= 


P-a^-R^D 

4E^I 


■(6  fi  -  8  sin(^)  +  sin(2  /i)) 


2.26.a 


PRb  (E-K) 
(Ea^  -  Kb^) 


■db  +  PR 


—  oldb 

\db  / 


2.26.b 


and 


Za(RPclP)  +  (1.5PDclR) 


2.26.C 


According  to  the  principle  of  virtual  work  defined  in  Ref  11,  the  change  in 
external  work  is  equal  to  the  change  in  elastic  strain  energy  {dW  -  dU)  if  a  configuration 
satisfies  geometric  constraints  and  remains  continuous.  Assuming  that  the  Bourdon 
satisfies  these  conditions. 


+  dU^  2.27 

where  the  incremental  external  work  acting  on  the  tube  is  defined  as 

dW*PdV  2.28 

or 
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2.29 


(P  +  dP)dV=dU  +dU. 

^  ^  xs  b 

Making  the  appropriate  substitutions  into  equation  2.29,  the  resulting  equation  is  written 
as 


(a^  +  b^)  E-2K  b^  „ 

(P  +  dP)-7ca-^ - i - db= 

Ea^-Kb^ 


2  a  E^  I  C  db^ 
1 

1  -  V 


rK 


-0 


^C^  dn'‘^  u  +  C^  dn''*  u+ C^  dn'^  u^  du  + 


Pa^RD 

4E'I 


•(l-v'>(6-p 


8  sir^  +  sir(2  p))  x 


PRb-(E-K) 

(Ea^-Kb^) 


/  d  ^  ^ 

Id  \  3PD 

a— D  +  D 

db+p.R.  db  +  PRdP+  dR 

\  da  1 

\db  /  ^  , 

2.30 


which  is  quadratic  in  db  where  dR  is  also  unknown.  See  Appendix  B  for  derivation 
details. 


4.  Tip  Travel 

The  equation  for  tip  travel  was  derived  by  making  the  assumption  that  the  length 
of  the  Bourdon  remains  constant  for  any  given  pressure  and  by  applying  Castigliano’s 
second  theorem.  The  assumption  for  constant  length,  L,  written  as 
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^iR=L 


2.31 


implies  that 


pdR  +  Rdp=0 


2.32 


or  that  incremental  tip  travel  dR  can  be  calculated  using  the  expression 


dR= 


-Rdp 


2.33 


Whereas  R  and  //  are  known  quantities,  the  variable  d/i  is  undefined.  Referring  to 
Figure  6, 

2.34 


by  geometry  where  ^  ,  the  angle  of  rotation  of  the  tip  due  to  the  applied  dummy  moment, 
is  described  by  the  equation 


rL 


M 

E^-I 


-'0 


—  M  dx 
^dM. 


2.35 


which  is  the  application  of  Castigliano’s  second  theorem  [Ref  11].  When  the  dummy 
moment,  M  is  applied,  the  equation  for  the  net  moment  becomes 

— 5 

Ms-P  a  R-D  (l  -  cos;S)  k-  M.  k  2.36 
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Figure  6.  Bourdon  length  segment  and  tip  deflection  for  the  application  of 
Castigliano’s  second  theorem 


Substituting  equation  2.36  into  equation  2.35  and  setting  the  dummy  moment  equal  to 
zero,  the  expression  for  tip  rotation  is  obtained  by  integrating 


,_PaRD 


(fi  -  sinpi) 


2.37 


In  this  formulation,  tip  displacement  is  calculated  iteratively.  For  this  model,  knowing  ^ 
defines  which  then  allows  the  determination  of  incremental  tip  travel  dR.  Once  dR  is 
known,  db  is  found  by  solving  equation  2.29  and  da  is  calculated  using  equation  2.8. 
Subsequent  values  of  dR  are  then  calculated  using  incremented  values  of  a  and  b. 


D.  SIMULATION  PROGRAM  AND  ANALYTICAL  RESULTS 

1.  Source  Code  Description  and  Limitations 

Considerable  effort  was  expended  in  expressing  all  of  the  equations  in  Chapter  II 
as  functions  of  a,  b,  k,  E,  and  K  to  simplify  the  programming  portion  of  this  research 
effort.  Having  derived  the  requisite  equations  to  determine  tip  travel  and  sectional 
deformation,  a  program  was  written  in  FORTRAN  to  iteratively  calculate  solutions  for 
user-defined  Bourdon  properties,  geometries,  and  operating  pressures.  The  source  code, 
listed  in  Appendix  E,  defines  all  program  variables,  prompts  the  user  for  input,  and  then 
calls  seven  subroutines  to  iteratively  calculate  tip  travel  and  corresponding  changes  in  tip 
geometry. 

The  topics  of  temperature  effects  and  the  standard  methods  of  compensation 
discussed  in  Chapter  I  were  not  addressed  in  this  attempt  to  develop  an  analytical 
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solution  for  Bourdon  tip  displacement.  The  analytical  model  presented  in  this  paper  only 
calculates  tip  travel  and  cross-sectional  deformation  as  functions  of  pressure. 


2.  Analytical  Results 

Execution  of  the  source  code  Program  Bourdon  was  very  sensitive  to  the 
geometric  relationships  between  a,  b,  and  R  and  worked  best  for  a:b  ratios  near  2;1  and 
for  0  <  <  34  mm  {Q<R<  1.34  in).  The  program  ceased  to  work  beyond  one  iteration 

once  R  exceeded  33  mm. 

The  following  table  indicates  analytical  results  for  fictitious  Bourdon  elements 
each  fabricated  from  phosphor  bronze  having  an  elastic  modulus  of  110  GPa  (16  x  10® 
psi)  and  a  Poisson’s  ratio  of  0.34; 


Table  1.  Analytical  Solutions  Generated  from  Program  Bourdon 


Bourdon 

Element 

Major  Axis 

(2a) 

[nim/in] 

i 

Minor  Axis 

(2b) 

[mm/in] 

i 

Radius  of 

Curvature 

(R) 

[rran/in] 

Wall 

Thickness 

(t)  [mm/in] 

Angle 

Describing 

Tube 

Length  (u) 

[deg/rad] 

Max. 

Applied 

Pressure 

(P) 

[MPa/psi] 

Tip  Travel 

(dR) 

[mm/in] 

1 

10.000  / 

5.000  /  1 

30.000  / 

250/ 

1.00  / 

1.54/ 

1 

1 

0.393 

0.197 

1.181 

i 

4.363 

145.00 

0.06 

2 

10.000  / 

5.000  / 

30.000  /  1 

0.100  / 

250/ 

1.00/ 

41.54  / 

0.393 

0.197 

1.181 

0.004 

4.363 

145.00 

1.64 

3 

10.0/ 

5.000  / 

30.000/ 

immiii 

330/ 

1.00/ 

1.38/ 

0.393 

0.197 

1.181 

0.012 

5.760 

1 

i 

145.00 

0.05 

4 

10.0/ 

5.000  / 

20.000  / 

0.300  / 

250/ 

1.00  / 

0.46  / 

0.393 

0.197 

0.787 

0.012 

4.363 

145.00 

0.02 

5 

10.0/ 

5.000  / 

50.000  / 

mm 

250/ 

1.00/ 

7.12/ 

0.393 

0.197 

1.97 

0.012 

4.363 

145.00 

0.28 

6 

14.000  / 

5.000  / 

30.000  / 

0.100  / 

250/ 

1.00/ 

112.41  / 

0.551 

0.197 

1.181 

0.004 

4.363 

145.00 

4.43 

7 

14.000  / 

5.000  / 

20.000  / 

250/ 

1.00/ 

33.31  / 

0.551 

0.197 

0.787 

4.363 

145.00 

1.31 

8 

7.500  / 

5.000  / 

20.000  / 

250/ 

1.00/ 

4.32/ 

0.295 

0.197 

0.787 

4.363 

145.00 

0.17 

9 

7.500  / 

5.000  / 

30.000  / 

250/ 

1.00/ 

14.57  / 

0.295 

0.197 

1.181 

4.363 

145.00 

0.57 

Looking  at  the  results  of  the  analytical  model,  some  of  the  tip  deflections  are 
excessively  large  while  others  are  of  more  reasonable  magnitudes.  These  extremes  in 
output  are  currently  unexplained  and  under  investigation.  While  the  validity  of  the 
analytical  model’s  results  is  unknown,  an  important  observation  is  made  concerning  the 
trends  of  the  output.  As  stated  previously  in  Chapter  I,  Bourdon  elements  having  flatter 
cross-sections  and  thirmer  walls  exhibit  larger  tip  displacements  than  rounder,  thicker 
elements  for  the  same  applied  pressure.  The  results  presented  in  Table  1  show  trends 
that  agree  with  Bourdon  displacement  theory. 

Comparing  Bourdon  elements  1  and  2  from  Table  1,  element  2  has  a  thinner  wall 
than  element  1  and  consequently  experiences  significantly  more  tip  travel.  Looking  at 
the  major  axes  of  elements  2,  6,  and  9,  element  6  has  the  flattest  cross-section  followed, 
in  increasing  order,  by  elements  2  and  9,  respectively.  As  expected,  element  6  exhibits 
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the  greatest  tip  movement  and  element  9  exhibits  the  least  when  all  other  parameters  are 
held  constant  between  the  three  specified  elements.  Similar  comparisons  can  also  be 


made  between  other  elements  in  Table  1. 


III.  EXPERIMENTAL  AND  NUMERICAL  ANALYSES  OF  A 

BOURDON  TUBE 

A.  EXPERIMENTAL  ANALYSIS 

Because  analytical  studies  are  based  on  assumptions,  it  is  good  practice  to 
compare  analytical  results  with  experimental  data.  Doing  so  verifies  the  validity  of  the 
mathematical  model  and  its  underlying  assumptions.  For  the  model  developed  in  this 
thesis,  a  good  correlation  between  analytical  results  and  experimental  data  implies  the 
successful  application  of  this  analytical  model  to  future  designs  of  Bourdon  gauges. 

The  experimental  results  presented  in  Tables  2  and  3  were  obtained  by  two  naval 
facilities:  the  calibration  laboratory  at  the  Regional  Repair  Center  (RRC),  Norfolk,  and 
the  Mechanical  Instrument  Repair  Calibration  Shop  (MIRCS)  in  USS  SHENANDOAH 
(AD-44),  homeported  in  Norfolk.  The  requested  information  was  faxed  to  SIMA, 
Norfolk  and  was  then  disseminated  to  the  aforementioned  facilities.  Pertinent 
dimensions,  elastic  element  information,  and  Bourdon  tip  travel  results  from  the 
experimental  studies  are: 
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Table  2.  Experimental  Data  from  RRC 


Gauge 

Company 

and  Range 

[psi] 

Bourdon 

Material 

Major  Axis 

(2a) 

[in/mm] 

Minor  Axis 

(2b) 

[in/mm] 

Radius  of 

Curvature 

(R)  [in/mm] 

Applied 

Pressure  (P) 

[psi/MPa] 

Tip  Travel 

(dR) 

[in/mm] 

Weksler 

Alloy  Steel 

0.5737  / 

0.2220  / 

1.361  / 

300.00  / 

0.25/ 

0-300 

14.5720 

5.6388 

34.5694 

2.07 

6.35 

Mensor 

Stainless 

0.5559  / 

0.2936  / 

3.0843  / 

1000.00  / 

0.4375  / 

0-1000 

Steel 

14.1199 

7.4574 

78.3412 

6.89 

11.1125 

Heise 

Stainless 

0.4108  / 

0.2758  / 

3.5129  / 

3000.00  / 

0.25/ 

0-3000 

Steel 

10.4343 

7.0053 

89.2277 

20.68 

6.35 

Table  3.  Experimental  Data  from  MIRCS 


Bourdon 

Material 

Major  Axis  (2a) 

[in/mm] 

Minor  Axis  (2b) 

[in/mm] 

Radius  of 

Curvature  ( R ) 

[in/mm] 

Applied 

Pressure  (P) 

[psi/MPa] 

Tip  Travel 

(dR) 

[in/mm] 

Stainless 

0.5/ 

0.115/ 

0.771  / 

500/ 

0.131  / 

Steel 

12.7 

2.921 

19.571 

3.45 

3.327 

Beryllium 

0.5/ 

0.190  / 

1.033  / 

600/ 

0.129  / 

Copper 

12.7 

4.826 

26.236 

4.14 

3.277 

One  key  parameter  missing  from  Tables  2  and  3  is  the  Bourdon  wall  thickness  of 
each  gauge.  Whereas  the  major  and  minor  diameters  were  easily  measured,  wall 
thicknesses  were  neither  measurable  nor  available  from  the  respective  manufacturers  due 
to  the  proprietary  nature  of  the  information. 
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B.  FINITE  ELEMENT  ANALYSIS 


One  rapidly  growing  approach  to  engineering  problem  solving  is  the  numerical 
technique  of  finite  element  analysis.  These  techniques  (1)  discretize  a  domain  with  a 
finite  number  of  finite  elements  and  a  finite  number  of  nodes  per  element,  (2)  introduce 
nodal  points  and  appropriate  shape  functions  to  represent  the  trial  solution  in  each 
element,  and  (3)  transform  the  resulting  set  of  differential  equations  into  a  series  of 
matrix  equations  to  change  a  system  from  an  infinite  degree  of  freedom  system  to  a  finite 
degree  of  freedom  system  of  appropriate  size. 

Several  powerful  finite  element  software  packages  exist  today  that  compute 
solutions  for  complicated  engineering  problems.  For  the  purpose  of  this  research,  the 
ability  to  build,  constrain,  and  pressurize  a  Bourdon  tube  using  computer  graphics  greatly 
enhanced  the  study  of  Bourdon  behavior.  As  powerful  as  this  tool  can  be,  however,  this 
type  of  analysis  did  not  preclude  the  effort  to  conduct  laboratory  experiments  but  rather 
supplemented  the  study  of  Bourdon  behavior.  Once  the  models  were  constructed  and 
drawn  to  match  the  dimensions  and  material  properties  of  the  Bourdons  used  in  the 
experimental  studies,  tip  displacement  results  were  quickly  determined. 

The  finite  element  solutions  presented  in  this  paper  were  generated  using  I-DEAS 
Simulation  software,  software  that  builds  geometry  and  finite  element  models,  analyzes 
the  models,  and  evaluate  the  results.  Each  of  the  finite  element  models  presented  in  this 
paper  was  constructed  using  the  following  parameters:  (1)  major  axis  2a,  (2)  minor  axis 
2b,  (3)  wall  thickness  t,  (4)  radius  of  curvature  R,  (5)  Bourdon  arc  length  //,  (6)  elastic 
modulus  E,  (7)  Poisson’s  ratio  v,  and  (8)  applied  pressure  P. 
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Assumptions  for  wall  thickness  and  arc  length  were  made  based  on  design 
information  presented  in  Refs.  1  and  5.  Boundary  conditions  were  applied  at  one  end  of 
each  Bourdon  model.  The  socket  or  pressure  end  was  clamped  to  prevent  rotation  and 
translation  in  the  x,  y,  and  z  directions.  The  tip  end  of  each  model  was  sealed  but 
otherwise  left  unconstrained  to  allow  maximum  freedom  of  movement.  Internal 
pressures  corresponding  to  experimental  test  pressures  were  applied  to  the  walls  and  end 
cap  (tip)  of  each  finite  element  model. 

The  numerical  models  were  analyzed  using  linear,  thin-shelled  elements  whose 
thicknesses  were  modified  to  match  the  estimated  wall  thicknesses  of  the  Bourdon 
elements  listed  in  Table  1.  The  quality  of  the  mesh  for  each  element  was  established 
using  criteria  for  distortion.  On  a  scale  of -1.0  to  1 .0,  an  element  is  perfectly  square  for  a 
distortion  value  of  1.0.  Meshes  exhibiting  distortion  values  less  than  zero  are  very  bad 
and  should  be  avoided  whereas  distortion  values  greater  than  0.4  are  acceptable.  In 
modeling  the  Bourdon  elements  using  I-DEAS,  all  of  the  meshes  exhibited  distortion 
levels  greater  than  0.4.  The  elements  that  were  distorted,  however,  consistently  occurred 
along  the  length  of  the  element  at  the  major  axis. 

The  following  table  summarizes  numerical  solutions  for  tip  displacement  in 
Bourdon  elements  having  the  same  geometries  and  material  properties  as  those  presented 
in  Chapter  1 1, Table  1; 
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Table  4.  Numerical  Solutions  for  Table  1  Bourdons 


Bourdon 

Element 

Major  Axis 

(2a) 

[mm/in] 

Minor  Axis 

(2b) 

[mm/in] 

Radius  of 

Curvature 

(R) 

[mm/in] 

Wall 

Thickness 

(t)  [mm/in] 

Angle 

Describing 

Tube 

Length  (g) 

[deg/rad] 

Max. 

Applied 

Pressure 

(P) 

[MPa/psi] 

Tip  Travel 

(dR) 

[mm/in] 

1 

10.000/ 

5.000  / 

30.000  / 

250/ 

1.00/ 

3.26/ 

0.393 

0.197 

1.181 

4.363 

145.00 

0.13 

7 

14.000  / 

5.000  / 

20.000  / 

250/ 

1.00/ 

42.8/ 

0.551 

0.197 

0.787 

4.363 

145.00 

1.69 

8 

7.500  / 

5.000/ 

20.000  / 

0.100  / 

250/ 

1.00/ 

8.50/ 

0.295 

0.197 

0.787 

0.004 

4.363 

145.00 

0.33 

Bourdon  elements  1,  7,  and  8  were  chosen  from  Table  1  at  random  for  modeling 
in  I-DEAS.  As  for  the  case  of  the  analytical  model,  the  validity  of  the  magnitudes  of  tip 
deflection  generated  by  the  finite  element  technique  is  also  uncertain  without 
comparative  data  from  gauge  manufacturers  or  more  extensive  experimental  work. 

Comparing  the  results  from  Tables  1  and  4,  the  numerical  solutions  follow  the 
same  trends  as  those  from  the  analytical  model  but  at  greater  magnitudes  of  tip 
deflection.  Tip  travel  is  53%,  22%,  and  49%  greater  for  I-DEAS  elements  1,  7,  and  8, 
respectively,  than  for  the  same  analytical  elements  in  Table  1. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  objective  to  analytically  derive  a  model  to  predict  Bourdon  behavior  was  a 
very  challenging  one.  The  very  mechanism  by  which  Bourdons  move,  as  described  in 
Chapter  I,  indicates  the  level  of  complexity  of  this  engineering  problem.  The  simple 
beam  and  thin-shell  theories  and  equations  used  in  the  development  of  this  analytical 
model  for  tip  travel  in  a  Bourdon  tube  have  produced  solutions  whose  validity  is  yet  to  be 
determined. 

The  task  of  comparing  analytical,  numerical,  and  experimental  results  was 
hampered  by  insufficient  experimental  data,  questionable  numerical  solutions,  and  forced 
assumptions  concerning  Bourdon  geometry  and  material.  The  fact  that  Program 
Bourdon  produced  output  that  followed  expected  trends  is  encouraging.  The  fact  that  the 
code  runs  for  more  than  one  iteration  only  for  certain  Bourdon  geometries  warrants 
continued  investigation. 

The  experimental  data  collected  from  RRC  and  MIRCS  was  very  limited  due,  in 
some  part,  to  the  underspecification  of  the  requested  information.  Ideally,  the 
dimensions  of  the  elastic  elements  should  have  come  from  the  gauge  manufacturers,  but 
as  stated  previously,  that  information  is  closely  guarded  by  the  companies  that  design  and 
manufacture  the  gauges.  In  addition,  the  tip  travel  for  each  bourdon  element  should  have 
been  measured  over  the  operating  range  of  the  gauge  and  not  just  at  the  maximum 
allowable  pressure. 
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The  accuracy  of  the  numerical  results  generated  in  I-DEAS  was  questionable  due 
to  the  excessive  magnitudes  of  tip  deflection  for  the  applied  pressures.  Because  the 
experimental  and  numerical  results  for  tip  travel  were  each  somewhat  inconclusive  and 
because  the  analytical  results  were  constrained  by  the  geometric  relationships  of 
dimensions  a,  b,  and  R,  the  comparative  study  of  results  yielded  little  useful  information 
for  us  to  be  able  to  state  with  conviction  the  accuracy  of  the  analytical  results. 
Consequently,  the  validity  of  the  analytical  model  remains  under  investigation. 


B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

The  research  results  presented  in  this  paper  are  rather  limited  implying  that  a 
considerable  amount  of  work  still  needs  to  be  done  to  develop  an  acceptable  analytical 
model  for  Bourdon  behavior.  The  following  suggestions  are  made  to  continue  the 
research  effort: 

1.  Peruse  Chapter  VII  (Tubular  Manometric  Springs)  of  Ref  12  for  its  detailed 
discussion  of  Bourdon  tubes  and  compare  the  assumptions  and  derivations  presented  in 
this  thesis  to  the  theories  presented  in  Elastic  Elements  of  Instruments.  This  book  is 
highly  recommended  by  the  Heise  Division  of  Dresser  Instruments  as  an  excellent 
resource  in  developing  an  analytical  model  for  Bourdon  behavior. 

2.  As  a  means  of  eliminating  the  geometric  constraints  of  the  analytical  model, 
modify  the  model  using  less  restrictive  geometric  assumptions  by  independently 
considering  variable  Bourdon  length,  circumference,  and  angle  f 
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3.  Collect  more  thorough  experimental  data  covering  a  much  broader  operational 
range  of  pressures  to  better  analyze  and  predict  Bourdon  behavior  and  eliminate  as  many 
assumptions  as  possible  concerning  Bourdon  geometry  by  obtaining  comprehensive 
measurements. 

4.  Enhance  the  analytical  model  to  include  a  stress  analysis  to  predict  failure 
from  fatigue  due  to  cyclic  loading  or  from  overpressurization  and  exploit  the  ability  of 
finite  element  software  packages  to  calculate  stresses  and  strains  within  finite  element 
models  for  comparison  purposes. 
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APPENDIX  A.  GEOMETRIC  RELATIONSHIPS 


Geometric  relationships  are  derived  in  great  detail  in  this  appendix.  These 
relationships  ultimately  serve  as  the  building  blocks  for  the  formulation  of  the  analytical 
model  for  Bourdon  behavior.  Derivations  for  da,  the  change  in  the  semi-major  axis,  dV, 
the  volumetric  change  of  the  tube  of  unit  length,  and  d(l/p),  the  change  in  the  curvature 
of  the  cross-section,  are  presented  in  this  appendix. 

Referring  to  Figure  4,  the  equation  for  the  elliptical  cross-section  is 


and  the  area.  A,  is 

AaTrab 

The  volume  of  the  elliptical  tube  with  length  unity  is 
VssTr-ab 

and  the  change  in  volume  is 
<JV5S5r(adb  +  bda) 

The  circumference,  1,  of  an  ellipse  is  defined  as 


A.2 


A.3 


A.4 
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where  ds  is  an  incremental  length  of  the  circumference  defined  by  the  radius  of  curvature 
of  the  cross-section  p  and  the  angle  (f>.  Referencing  Figure  4  for  small  d^. 


ds 


4 


dx^  +  dz^ 


A.6 


For 


x=bcos(^ 


A.7 


and 


z=  a-  sini/> 


A.  8 


the  following  expressions  are  true: 

dx^ =b^-  sin^-  d^^  A.9 

and 

dz^=a^cos^<^d<)!»^  A.  10 


Substituting  equations  A.9  and  A.  10  into  equation  A.6  and  dividing  both  sides  by  the 
square  root  of 


ds=a- 


A.11 


Substituting  in  the  trigonometric  identity  cos^(f>  =  1  -  sin^t/i  and  rearranging  terms. 


sin  +  l  d(^ 


A.  12 


The  new  variable  ^is  defined  by  the  expression 


2 

a 


A.  13 


on  the  interval  0  <  A:  <  1.  The  incremental  arc  length  ds  is  now  written  as 


ds*a' 


-  k^-  sin*-  </>•  d<^ 


A.  14 


so  that  the  circumference  becomes 


X 

~2 


A.  15 


and  then 

1=4- a- E 


A 16 
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where  E  is  the  complete  elliptic  integral  of  the  second  kind  defined  by  equation  2.2.  For 
the  assumption  that  /  is  constant,  dl  =  0.  Taking  the  derivative  of  equation  A.  16  and 
equating  it  to  zero  yields  the  following  expression: 


a 


(in' 

•dk 

\dk  i 

+  E-  da®0 


A.  17 


Solving  for  dk,  the  derivative  of  both  sides  of  equation  A.  14  is  taken  and  written  as 


kdk- 


-j  (a  db  -  b  da) 
& 


A.18 


then  as 


b 

dk= - (b  da  -  a  db) 

ka^ 


A.  19 


and  eventually  as 


dk= 


b  da  -  a  b  db 


2  /  2 
•,^a  - 


A.20 


after  making  the  substitution  for  k  defined  by  equation  A.  13.  From  Ref  9, 


d  —  E  -  K. 
_E«= - 

dk  ^ 


A.21 
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where  K  is  the  complete  elliptic  integral  of  the  first  kind  defined  by  equation  2.1. 
Substituting  equations  A.20  and  A.21  into  equation  A.  17  yields 


-m 


/b  da  -  a  db\ 
\  k-a^  / 


Eda=o 


A.22 


which  simplifies  to 


(E  b^  -  K  b^)  da  +  (K  a  b  -  E  a  b)  db  +  (E  a^  -  E  b^)  da 

a  -  b 


A.23 


and  eventually  to 


ab(E-K) 

E-a^-K-b^ 


•db 


A.24 


the  ultimate  expression  of  the  change  in  the  semi-major  axis  written  as  a  function  of  the 
change  in  the  semi-minor  axis.  Substituting  this  expression  into  equation  A.4,  the 
volumetric  change  in  the  Bourdon  may  be  written  as  a  fimction  of  db  as  follows: 


dV=7C- 


a-  db  +  b- 


ab(E-K) 
E  a  -  K  b^ 


After  some  algebraic  manipulation,  the  change  in  volume  is  written  as 


A.25 


dV=!?ra- 


(a^+b^)  E-2  K  b^  „ 
- db 

Ea^-K-b^ 


A.26 
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Referencing  Figure  4,  the  radius  of  curvature  p  of  the  cross-section  is  defined  by  the 


relation 


and  the  curvature  is  defined  by  the  expression 


^  /dz' 


The  objective  now  is  to  write  curvature  as  a  function  of  the  parameters  a,  h,  and  (f). 
Recalling  equation  A.1,  its  derivative  is 


X-  dx  z-  dz 


Solving  for  dz/dx. 


dz  -a^  x 
dx  b^.2 


so  that 


dz\^  a^  x^ 


The  next  objective  is  to  write  dz^/db^  in  terms  of  a,  b,  and  (f).  By  definition, 


With  the  substitution  of  equation  A.30,  equation  A.32  becomes 


J  ^  2  \ 

d  z  d 

dx^  A.3: 


or 


dx^  dx  \zj  A.  34 


which,  when  the  derivative  is  taken,  becomes 


—-'A 


A.35 


Substituting  equation  A.30  into  equation  A.35  and  simplifying,  equation  A.35  may  be 
rewritten  as 


d^  z 

dx^  b^z 


• 

1  2  \ 

^  X 

z 

f-a  x 

^b^zj 

A.36 


and  eventually  as 


Now  substituting  equations  A.30  and  A.37  into  equation  A.28,  curvature  is  expressed  as 


1 

P 


-a 


3 

b  Z 


1 


4  2\  1-5 

a  x  \ 
b  z  / 


A.38 


which,  after  a  few  algebraic  manipulations,  becomes 


1 

P 


■a  b 


(b  z  +  a  x  ) 


.2,1.5 


A.39 


Substituting  equations  A.9  and  A.  10  into  equation  A.39,  curvature  is  written  as 


1  -a  b 

jO  /■U+  2  2  Jl  ♦  u2  2  J 

^  (b  a  Sin  9  +  a  b  COS  9) 


A.40 


or  as 


1  -a  b 

jO  /i_2  -Si  2  2  1 A 

^  (b  Sin  '9  +  a  cos  f) 


A.41 


Recalling  the  trigonometric  identity  cosr^  =  J  -  sin^ij),  and  dividing  the  numerator  and 
denominator  by  , 


K  fb‘ 


■  sin^-  <j(>  +  (1  -  sin^-  <!>)  ■  ^ 


and  ultimately 


Considering  the  absolute  value  of  the  curvature 


^  a^  (l  -  k^  sin^'^)^ 


which  now  correlates  to  the  definition  of  curvature  for  an  elliptical  cross-section  as 
defined  in  Ref  13.  The  change  in  curvature  may  now  be  written  as 


where 


P  a*  (l  -  k^  sin^  i^)^ 
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\ 

—  pl-db  + 
,db  / 


A.47 


The  derivative  is  taken  with  respect  to  k  because  A:  is  a  function  of  the  dimensions  a  and  b 
as  defined  by  equation  A.  13.  The  partial  derivatives  of  p  with  respect  to  a,  b,  and  k  are 


-  k^  sin^-(^)^'^  da 


A.47.a 


I  —  jO  I  •  db= — •  ( 1  -  k^-  sin^-  <j))  ■  db 
\db  /  b^ 


A.47.b 


■dk* 


-3a*k-  dksin^- 
b 


A.47.C 


so  that  equation  A.47  becomes 


—■  dbj  (1  -  k^-  sin^-  <p) 


a^k-  dksin^- 


<!> 


Jl  -  k^  sin^  i][> 


Recalling  equations  A.  13  and  A.  19,  the  expression  kdk  may  be  written  as 


A.48 


k-  dk= —  (b-  da  -  a-  db) 
a^ 


A.49 
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Making  the  substitutions  for  da  and  kdk  and  simplifying,  equation  A.48  becomes 


,  a"  (2b^-a")  E-Kb'  ,,  ,,  .  j  ^ 

dpss — •  -  (1  -  k  sin  -  9)  db 

Ea^  -  Kb^ 


3  E  sin^  <^- ^ -  Jl  -  k^- sin^- ij[>- db 

lEa^  -  Kb^J 


which  is  even  further  simplified  by  manipulating  the  following  terms; 


2  .-2  j_  t  1  ,2  .2  1  1  2  _■  2 


sin  -  k  sin  - -  k  sin  (^  ((l  -  k  sin  1^)  -  l) 


Ja  L  1 2  •  2  1 2  •  2  ,  ,2-2 


sm  ^  -  k  sin  -^ss — [(l  -  k  sin  •^)  -  -  k  sin  ip 


Making  the  substitution  of  equations  A.  13  and  A.  52  into  equation  A.  50  ,  the  expression 
for  the  change  in  the  radius  of  curvature  becomes 


3‘ a  -E  5  T 

4o= - -^1  -  k^- sin^- ijl>- db  + 


Ea^  -  Kb^ 


(2  b“  -  a')  E  -  K  b"  3  a"  E  „  ,  i  ,  ,,i.5  ^ 

- (1-ksin®)  db 

b^  Ea^  -  Kb^  Ea"  -  Kb^ 


which,  with  a  little  more  manipulation,  becomes 
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Ea^  -  Kh^ 


Now  substituting  equations  A.46  and  A.54  into  equation  A.45,  the  change  in  curvature  is 


written  as 


^  Ca%b^)E  2  2 

df-  = - - - 3E-  - ~  +  K  (l-k^  sin^- 

P!  /(l-k^sm^-.^f'CEa^-Kb^)l  b^ 
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APPENDIX  B.  INCREMENTAL  STRAIN  ENERGY  DUE  TO 
CROSS-SECTIONAL  DEFORMATION 


In  general,  elastic  strain  energy  for  a  beam  subjected  to  transverse  loading  is 


expressed  as 


So  as  not  confuse  the  symbols  for  elastic  modulus  and  elliptic  integrals  of  the  second 
kind,  elastic  modulus  for  a  straight  beam  will  be  designated  E‘ ,  elastic  modulus  derived 
from  the  theory  of  thin  shells  will  be  designated  E,  ,  and  elliptic  integrals  of  the  second 
kind  will  be  designated  E  for  all  remaining  calculations.  For  an  initially  curved  beam, 
the  moment  Mis  defined  by  the  expression 


M=E  -I- 


Incremental  strain  energy  in  the  bourdon  due  to  a  change  in  the  cross-section  is  written  as 


dU  =4-  - ds 

2E-I 

1 


where  Mis  now  defined  as 
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B.4 


M=EId 

1 


Making  the  appropriate  substitutions  for  MandE/,  equation  B.3  becomes 


dlJ 


2-FI 


XS  2 

1  -  V 


% 

1 


-'0 


where,  for  a  unit  length,  the  area  moment  of  inertia  is  defined  as 


B.5 


12 


B.6 


and  t  is  the  wall  thickness  of  the  Bourdon  element.  The  integral  term  is  found  by 
squaring  equation  A.45  which  results  in  the  expression 


2 


3E- 


(a^  +  b^)  E 
b^ 


+  K 


(1 


k^-  sin^ 


db 


B.7 


which  upon  squaring  becomes 
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Rewriting  equation  B.5  in  terms  of  the  constants  Cy  through  C4  and  substituting 
equation  A.  14  for  ds,  the  equation  for  incremental  strain  energy  due  to  a  change  in  the 
cross-section  becomes 


dU 


xs 


2E^I 


1  -  V 


C  db  ■ 
3  1 


(1  -  k^  sin^-^)  (1  -  k.^  sin^  (^)  (l-k  sin  -«;6) 


■ajl 


k^-  sin^-  d<^ 


B.IO 

With  further  simplification,  incremental  strain  energy  due  to  sectional  deformation  is 
written  as 


db> 

1  -  V 


C 


-'0 


(1  -  k^  sin^-^)^^  (1  -  (1  - 


B.ll 

All  of  the  terms  in  the  denominator  of  the  integral  in  equation  B.  1 1  may  be  written  in 
terms  of  elliptic  fimctions  by  letting 


sin(«^)=ssn-(u) 


B.12 


and 
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cos(<j(>)=cn(u) 


B.13 


so  that 

d(  sm<^) = cosi^-  dj>st  d(  sn(u) ) = cn(u)  •  dtt(u)  ■  du 
Solving  for  dtj), 

d<^=dn(u)du 
In  addition, 

1  -  k^  sin^  (<^)=l  -  k^  sn^  (u) 

where 

1  -  k^  sn^-(u)  =  dn^  Cu) 

by  definition  in  equation  2.4. b.  Equation  B.ll  is  further  simplified  by  evaluating  the 
integral  in  parts  and  by  making  substitutions  with  equations  B.15  and  B.17.  The  first 
term  in  the  integral  is  written  as 

■  % 

1  rK 

( 1  -  k^-  sin^-  d<^a  ( dn^-  (u) )  •  dn-  (u)  du 

'0  -0 


B.14 


B.16 
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(<in'^  (u)  dtt  (u))  du=  dn''^  (u)  du 


Applying  equation  5. 13 1 .3  from  Ref  14  for  m  =  -4, 


dn'*  (u)  du= 


3(kV 


rK 

■K 

k^'  dn'  (u)  sn(u)  •  cn(u)  -  2-  (2  -  k^)  ■ 

2 

dn'  •  (u)  du  +  u  du 

-0 

'0 

Using  the  elliptic  function  identities  defined  in  Chapter  I  and  the  relationships 


(dn^  (u))  du=E 


as  defined  in  Refs  9  and  10,  equation  B.  18  is  rewritten  in  simpler  form  as 


■K  ,  fK 

dn‘*-(u)  du= - ^ dn  ^  (u)  du 

0  3(1 -k^)  3(l-k^).0 


Again,  applying  equation  5. 13 1.3  from  Ref  14  to  solve  for  the  final  integral  where 
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m  =  -2,  equation  B.22  is  further  simplified  to  the  expression 


■K 

Hn  ( 11 '^1"* 

'2E(2-k") 

-0  3(1 -k^) 

(l-k^) 

B.23 


where,  for  future  reference, 


-2  E 

dn  (u)  du= - 

-0  d-k') 


B.24 


Still  in  the  process  of  simplifying  the  equation  for  the  incremental  strain  energy  in 
Bourdons  due  to  changes  in  cross-section,  the  next  integral  to  be  evaluated  is 


7C 

T 


rK 


(1  -  k^  sin^  <;t>) 


(dn^  (u))  ^’^■dn  (u)  du 


-'0 


rK 


rK 


(dn'  (u)  dn  (u))  du= 


-'0 


dn'*  (u)  du 


Applying  equation  5.131.3  from  Ref  14  for  m  =  -6, 


B.25 
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■K 

dn''*  (u)  du= 

-0 


1 


k^  dn'^  (u)  sn(u)  cn(u)  -  4- (2  -  k^)- 

•K 

dn"^  (u)  du  +  3- 

■K 

dn'^-(u)  du 

- 

0 

0 

B.26 

Substituting  equations  B.20,  B.23,  and  B.24  and  applying  the  elliptic  function  identities 
defined  in  Chapter  I, 


•K 

'0 


dn'*  (u)  clu= 


4(2-k^) 
15- (1  -  kV 


2-E(2-k^) 
(1  -  k") 


3-E 

5(1 -kV 


B.27 


And  finally,  applying  the  same  procedure  to  the  denominator  of  equation  B.ll 
corresponding  to  the  constant  C4 , 


■  % 

7 

0 


rK 


(1  -  k^  sin^M^) 


( dn^-  (u) )  ^  dn-  (u)-  du 


K 


rK 


(dn'  •(u)  dn  (u))  du=! 


dn'*  (u)  du 


where,  applying  equation  5. 13 1 .3  from  Ref  14  for  m  =  -8, 


B.28 
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•K 

dn'®  (u)  dus= 

'0 


•K 

■K 

k^  dn‘^  (u)  sn(u)  cnCu)  -6  (2-  k^)- 

dn'  **•  (u)  du  +  5- 

dn'*  (u)  du 

- 

0 

0 

B.29 

Substituting  equations  B.20,  B.23,  and  B.27  and  applying  the  elliptic  function  identities 
defined  in  Chapter  I, 


■K 

0  2 

2E(2-k^) 

dn'®  (u)  du« 

24  (2 -k  )  5 

- 

(1-k^) 

-0 

105- (1  -  k^)  21- (1  -  k^) 

18E-(2-k^) 

35- (1  - 


B.30 


The  equation  for  incremental  strain  energy  due  to  changes  in  cross-section  may  now  be 
written  as 


dU  = 

xs 


2  a  E^  I  C  db^ 
1 

(1  - 


K 


-0 


dn'^  (u)  +C^  dn"®  (u)  +  dn’®-(u) 


du 


B.31 


where  the  integrals  of  drt'  (u),  dn^(u),  and  dn^  (u)  are  defined  by  equations  2.5,  2,6,  and 
2.7,  respectively. 
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APPENDIX  C.  INCREMENTAL  STRAIN  ENERGY 
DUE  TO  BENDING 


As  defined  in  Appendix  B,  the  total  strain  energy  due  to  bending  in  a  beam  is  defined  as 


2E'I 


where  M  represents  the  net  moment  acting  on  the  bourdon  as  a  result  of  pressurization 
and  E'  is  Young’s  modulus.  The  detailed  derivation  of  the  net  moment  acting  on  the 
elastic  element  represents  a  very  large  segment  of  this  appendix. 

Referencing  Figure  5,  the  elemental  area  of  the  bourdon  tube  is  defined  by  the 


expression 


dA=!(R  +  pcos(t))-d0ds 


Recalling  equation  A.  14,  equation  C.2  may  also  be  written  as 


dA=a  (R  +  pcos(t))  A|i  -  k  sin  ^  dO  dcj) 


The  force  acting  on  the  elemental  area  is  defined  as 


dF=PdA 


where 


Ps=P((cosQ!Cos6)i  +  (cosQ!sin0)j  +  since- k) 
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The  angle  a  describes  the  direction  of  the  normal  component  of  pressure  acting  on  the 
wall  of  the  tube  as  shown  in  Figure  5.  Substituting  equations  C.3  and  C.5  into  equation 
C.4,  the  elemental  force  acting  on  the  elemental  area  of  the  bourdon  is 

— 5  r 

I  2  2 

dF=P-((coso!  cos^)  i  +  (cosa  sin^)  ]  +  sinG!-k)  [a  (R  +  /O  cosi/))>]l  -  k  sin 

C.6 

where  the  corresponding  moment  arm  to  the  elemental  area  is  defined  as 

— J 

dr=-R(cos;Si  +  sinj5j)  +  (R  +  pcos<^)(cos^i  +  sin^j)  +/3sin^k 
Relating  the  angle  a  to  the  angle  ^  in  Figure  6, 

COS(t) 

COSa = — r==zz 

Ji  -  k^-sirf-<t> 
and 

b-  sinb 

SlltXS! — - 

a  ji  -  k^  sirf  cj) 

As  a  result,  equation  C.6  may  be  expressed  as 

— ^  b 

dFssP  a  (R  +  pcoscj))-  (cos<t)  cos0)  i  +  (cos(j)  siii0)  j  +  — sin|)k  •d(|)  d0 

L  ^  J 

CIO 
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The  moment  acting  on  the  elemental  area  due  to  an  elemental  force  is  defined  by  the 
equation 


dM»dr  X  dF 


C.ll 


Temporarily  factoring  out  Pa(R+ pcos(l))d<f)d 6  equation  C  IO  for  simplicity,  the 

vectorial  components  of  the  resulting  cross  product  of  equation  C.l  1  are: 


cos<j)  sii6  p  sint)  — •sin|)  (R  sir6  -  R  sirp  +  p  cos<t)  sir6) 
a 


■1 


C.ll.a 


cos<|)  cos0  p  sin|) - sin|)  (R  cos0  -  R  cosP  +  p  cos(j)  cos0) 

a 


•J 


C.ll.b 


and 

((Rcos0  -  RcosP  +  p  cos(|)  cos0)  cos<|)sii6  -  cos(t)cos0  (R  sii6  -  R  sii<3  +  p  cos«t)  sii6))  k 

C.ll.c 

For  the  sake  of  brevity,  the  cosine  function  will  be  abbreviated  by  “c”  and  the  sine 
function  will  be  abbreviated  by  “5”  from  this  point  to  the  end  of  Appendix  C.  Continued 
simplification  begets  the  following  equation  for  the  moment  acting  on  the  bourdon  tube: 


*2-7C 

*P 

r 

^  1)  b 

Pa(R+  pc<|))- 

pc<t)S(t)s0 - (R+  pc(|))s(t)s0  +  — Rs(l)sp 

L 

[a  a  J. 

» 0 

0 
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rp  ,  , 

b  b 

P  a  (R4- p  c(j))'  p  c<t)  S(t)  c6 - (R  +  p  c(t))  s(j)  c0  +  — R  s(l)  cp  jdGdtj) 

L  L  ^  ^  jJ 

*'  0  ^  0 


•2-71  1*3 

+  P  a  (R+ p  c(j))  R  c(|)  (c0  sP  -  s0  cp)  kd0d(j) 

•  0  -0 

C.12 

Factoring  back  in  the  expression  Pa(R+pcosip)d^0  and  evaluating  the  x-component  of 
the  equation  C.12  first. 


r2-7i  fp 

b  b 

P  a-  (R  +  p  c(t>)-  p  c(t)-s0 - (R+ p  c(|))  s0  +  — R-sP  -scj)  d0  d(t)*o 

[a  a  J 

•-0  » 0 

C.13 

because  sin^  is  an  odd  fimction  on  the  interval  0  to  27i.  Similarly  for  the  y-component  of 
equation  C.ll, 

r2-3r  fp 

b  b 

-P  a-  (R+p-c<t))-  p-c(t)-c0 - (R+ p-c(t))-c0  +  — RcP  -s(j)  d0  d(|)=o 

!_  a  a  J 

» 0  •'0 

C.14 

Evaluating  the  z-component  of  equation  C.ll, 
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PaR- 


(R  +  p  c(j))  c(t)  (c0  sp  -  s0  cp)  d0  d(j)= 


-'0  ►  0 


4Pa^R 

3 


■(i-cp> 


K-E 


+  (2E-K) 


C.15 

Therefore,  the  final  expression  for  the  moment  acting  on  the  walls  of  the  tube  is 


M= 


4Pa-R 


(i-cP)- 


K-E 


+  (2  E-  K) 


•k 


Again  referring  to  Figure  5,  the  net  force  acting  on  the  end  cap  of  the  bourdon  element 


F 

cap 


s-P-jta- 


b-j 


C 


and  the  corresponding  moment  arm  is 


dr=R(i  -  cp)  i-  Rsp  j 


C. 


Performing  the  cross  product. 


M  sdr  X  dF  =-P  R-7r-a  b  (1  -  k 


cap 


n  IQ 


Therefore,  the  net  moment  acting  on  the  tube  is  the  sum  of  the  moments  acting  on  the 
walls  and  end  cap  of  the  Bourdon  as  shown  in  next  equation: 


- j  _5  - j 

M  =M  +  M  =P  a  R  (l  -  c^e)- 

net  cap  '  ^ 


4  a 

K-E  ^ 

— —  +  (2  E-  K) 

-  TT  b 

3 

C.20 


For  simplicity,  the  new  variable  D  is  defined  as 


D=a- 


4  a 
3 


K-E 


+  (2-E-  K) 


Tih 

a 


C.21 


such  that 


- ^ 

M^^^=PaRD(i-c(3)k 


C.22 


Substituting  equation  C.22  into  equation  C.l,  the  total  bending  strain  energy  is  expressed 
as 


(P  a  R  D  (1  -  c^)y 
2E^-I 


•Rd^ 


'0 


C.23 


The  incremental  strain  energy  is  calculated  by  differentiating  Ub  with  respect  to  the 
dimensions  a  and  b,  the  applied  pressure  P,  and  the  radius  of  curvature  of  the  tube  R  such 
that 
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— U  IdR 
\dR  / 


C.24 


•db  + 


db  / 


•dP 


where 


,  P^aR^D 
■das: - 

4'E^-I 


■(6  (i  -  8'sin(/i)  +  sm(2-/i))- 


^a-— D+D 

\ 


■da 


C.25 


P^a^R^^D 

4-E^^I 


(6- pi  -  8^sin(/i)  +  sin(2 


C.26 


Pa^R^^D^ 

4E^^I 


■(6- pi  -  8  sm(/x)  +  sm(2 7i))  dP 


C.27 


and 


dR= - 

SE^^I 


■(6  pi  -  8^sin(^)  +  sin(2^/x))^  dR 


C.28 


Summing  equations  C.25  -  C.28  and  simplifying,  the  incremental  strain  energy  in  the 
bourdon  due  to  bending  may  be  expressed  as 
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Pa^R^-D 

4E^-I 


■{6  {i  -  8  sinC^)  +  sin(2-^)) 


X 


PRb(E-K) 

Ea^-Kb^ 


(alo  +  D 

\  da 


•db+PR 


+  RDdP  +  l.j(PDdR) 


C.29 


where 


iD=--  ^'  (^-E)  (3  (a^-b^)-2-a^) 
da  ^  (a^  -  b^)^ 


(2  E  -  K) 


C.30 


and 


2a^b(K-E) 
(a"  -  b')' 


-  % 


C.31 


as  a  result  of  taking  the  partial  derivatives  of  equation  C.21  with  respect  to  the  semi¬ 
major  and  minor  axes,  respectively. 
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APPENDIX  D.  DEVELOPING  EQUATIONS  FOR  TIP  TRAVEL 

Two  equations  are  used  to  develop  the  equation  for  Bourdon  tip  displacement. 
The  first  equation  is  based  on  the  assumption  made  in  this  analytical  study  that  the  length 
of  the  elastic  element  remains  constant  during  pressurization.  Based  on  this  assumption, 
the  following  is  true; 


/i-Rss  Constant  D.l 

and 

/X  dR  +  R  d^aO  D.2 


where  incremental  tip  travel,  dR,  is  expressed  as 


clR» 


-Rd^ 


D.3 


Referring  to  Figure  6,  the  variables  R  and  ju  are  measurable  quantities  but  d/n  must 
be  calculated.  One  way  to  do  so  is  to  apply  Castigliano’s  Second  Theorem  as  defined  in 
Ref  1 1  which  states  that 


il -/I-m]  dx 
I  ^dM.  j 


where M  is  an  imaginary  moment  and  is  Young’s  modulus.  This  equation  is  helpful 
because,  in  beam  theory,  9^  represents  the  change  in  the  slope  of  the  elastic  member  at 


the  point  of  an  applied  imaginary  moment.  Geometry  then  dictates  that 


thereby  establishing  the  second  equation  necessary  to  calculate  tip  deflection.  Modeling 
the  Bourdon  as  a  beam,  this  relationship  makes  sense  because  as  the  radius  of  curvature 
at  the  tip  of  the  Bourdon  becomes  larger  vvith  pressure,  the  change  in  slope  becomes 
larger  and  the  angle  n  describing  the  length  of  the  elastic  element  becomes  smaller. 
Therefore,  once  (f>i  is  known,  bourdon  tip  displacement  is  calculated  using  equation  D.3. 

Continuing  with  the  application  of  Castigliano’s  second  theorem,  a  dummy 
moment  M  is  applied  in  the  direction  chosen  positive  for  so  that  the  bending  moment 
in  equation  D.4  becomes 


M=s-P  a  R  D  (l  -  cosj6)-k-  M-  k 


D.6 


and 

d  , , 

- M=-l 

dM. 

1 

Making  the  appropriate  substitutions  into  equation  D.4,  setting  M,  equal  to  zero,  and 
integrating. 


,_PaE.D 


0 


(1  -  cosid)  dd= 


PaRD 

E^I 


■{}i  -  sinp;) 


D.8 
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A  positive  value  of  6!  means  the  tip  angle  of  rotation  has  the  same  sense  of  direction  as 
the  applied  dummy  moment.  Because  is  a  function  of  the  dimensions  a  and  b  (a 
function  of  b  because  Z)  is  a  function  of  b),  an  iterative  solution  for  tip  travel  incorporates 
the  cross-sectional  deformation  that  occurs  simultaneously  with  straightening. 
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APPENDIX  E.  FORTRAN  SOURCE  CODE 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

♦ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


Program  Bourdon 

This  program  calculates  tip  deflection  in  low  pressure,  C-type  Bourdon 
elements  for  user-defined  properties,  geometries,  and  pressures. 


Definition  of  variables: 

a  =  semi-major  axis  of  the  elliptical  thin-walled  tube  (ID) 
b  =  semi-minor  axis  of  the  elliptical  thin-walled  tube  (ID) 
t  =  wall  thickness  of  elliptical  tube 
R  =  tip  radius  of  curvature  of  the  elliptical  tube 
p  =  internal  pressure  applied  to  the  tube 

mixs  =  I,  the  area  moment  of  inertia  of  the  tube  for  cross-sectional 
deformation 

mib  =  I,  the  area  moment  of  inertia  of  the  tube  for  bending 

mu  =  arclength  of  the  elliptical  tube 

epr  =  elastic  modulus  of  the  tube  material 

nu  =  poisson's  ratio  of  the  tube  material 

rad  =  input  variable  for  the  initial  radius  of  curvature 

rin  =  initial  radius  of  curvature  expressed  in  meters 

da  =  incremental  change  in  a 

db  -  incremental  change  in  b 

dR  =  incremental  change  in  R 

dp  =  incremental  internal  pressure  applied  to  the  tube 

dmu  =  incremental  change  in  mu 

smk  =  k,  defined  as  the  square  root  of  l-(b/a)'^2 
bk  =  K,  Jacobian  Elliptic  Integral  of  the  1st  kind 
be  =  E,  Jacobian  Elliptic  Integral  of  the  2nd  kind 

cl  =  a  constant  in  the  expression  for  cape 
c2dn  =  the  product  of  a  constant  and  the  elliptic  integral  of 
dn(u)e-4  in  the  expression  for  cape 
c3dn  =  the  product  of  a  constant  and  the  elliptic  integral  of 
dn(u)e-6  in  the  expression  for  cape 

c4dn  =  the  product  of  a  constant  and  the  elliptic  integral  of 
dn(u)e-8  in  the  expression  for  cape 
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D  =  a  constant  in  the  net  moment  equation  expressed  as  a 
function  of  a,  b,  K,  and  E 
dDda  =  partial  derivative  of  D  wrt  a 
dDdb  =  partial  derivative  of  D  wrt  b 

capa  =  a  constant  expressed  as  a  function  of  parameters  a,  b, 

P,  and  R  used  in  Subroutine  Bndenrg 
capb  =  a  constant  expressed  as  a  function  of  parameters  a,  b, 

P,  and  R  used  in  Subroutine  Bndenrg 
cape  =  a  constant  expressed  as  a  function  of  dimensions  a  and 
b  used  in  Subroutine  Xsecenrg 

capj  =  an  expression  representative  of  incremental  external  work 
capl  =  the  ratio  of  capj  to  capa  used  in  the  expression  for  capy 
capx  =  the  ratio  of  cape  to  capa  used  in  the  final  expression 
for  db 

capy  =  an  expression  representing  the  difference  in  capb  and 
capl  used  in  the  final  expression  for  db 
capz  =  a  constant  expressed  as  a  function  of  parameters  a,  b, 

P,  and  R  used  in  Subroutine  Bndenrg 
thetai  =  bourdon  tip  rotation  due  to  applied  dummy  moment 

dV  =  volumetric  change  in  the  tube,  dV  =  vol*db 
duxs  =  change  in  internal  energy  due  to  a  change  in  the 
cross-section  of  the  elliptical  tube  as  a  result 
of  an  applied  internal  pressure 
dub  =  change  in  internal  energy  due  to  bending  of  the 
elliptical  tube 


implicit  real  *8  (a-h,o-z) 
real*8  mixs,  mib,  mu,  nu,  maj,  min 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 

+  capl,capx,capy 

open  (7,file='data',status- old') 

read  (7,*)maj  ,min,  wl,rad,mu,dp,pmax,epr,nu 

close  (7) 

Dimensions  converted  to  meters. 


a  =  maj/ 1000.0 
b  =  min/1000.0 


t  =  wl/1000.0 
rin  =  rad/1000.0 
R  =  rin 

call  capke 


pi  =  3.141592654 

vol  =  pi*a*(((a**2  +  b**2)*be  -  2.0*bk*b**2)/(be*a’''*2 
+  -  bk*b**2)) 

smk  =  sqrt(1.0  -  (b/a)**2) 

mi  =  (pi/4.0)*((a  +  t)*(b  + 1)**3  -  a*b**3) 


*  This  loop  increases  the  internal  pressure  to  a  pre-determined 

*  maximum  pressure  in  increments  of  dp  so  that  associated  changes 

*  in  the  variables  a,  b,  R,  and  mu  may  be  considered  in  the  problem. 

* 

Open(4,file='data',status='new') 

* 

Do  10  p=0.0,pmax  -  dp, dp 
call  capke 

* 

vol  =  pi*a*(((a**2  +  b**2)*be  -  2.0*bk*b**2)/(be*a**2 
+  -  bk*b**2)) 

♦ 

call  xsecenrg 
call  dercapd 
call  dummy 
call  tip 
call  bndenrg 

* 

capj  =  (p  +  dp)*vol 
capl  =  capj/capa 
capx  =  capc/capa 
capy  =  capb  -  capl 

* 

*  The  final  expression  for  db  is  written  as  a  quadratic  equation  for 

*  both  roots  of  the  equation. 

* 

db  =  (-capy  +  sqrt(capy**2  -  4.0*capx*capz))/(2.0*capx) 

*  db  =  (-capy  -  sqrt(capb**2  -  4.0*capx*capz))/(2.0*capx) 

da  =  (a*b*(be  -  bk)/(be*a**2  -  bk*b**2))*db 
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*  Variables  incremented  through  the  loop. 

* 

a  =  a  +  da 

b  =  b  +  db 

R  =  R  +  dR 

mu  =  mu  +  dmu 

smk  ==  sqrt(1.0  -  (b/a)**2) 

mixs  =  t**3/12.0 

mi  b=  (pi/4.0)*((a  +  t)*(b  + 1)**3  -  a*b**3) 

* 

10  continue 
* 

close  (4) 

* 

stop 

end 

* 

Subroutine  Capke 

* 

*  This  subroutine  calculates  K  and  E,  elliptic  integrals  of  the  first 

*  and  second  kind,  respectively.  K  and  E  are  functions  ofk  which  is 

*  itself  a  function  of  the  major  (2a)  and  minor  (2b)  axes  of  the 

*  ellipse. 

* 

implicit  real*8  (a-h,o-z) 
real*8  mixs,  mib,  mu,  nu 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj , 

+  capl,capx,capy 


ang  =  asin(smk) 
eps  =  l.Oe-8 
a0=1.0 
bO  =  cos(ang) 
cO  =  sin(ang) 
rmult=  1.0 
stor  =  c0**2 
bk=  1.5707963 
be  =  bk*(1.0  -  0.5*stor) 
89  if  (abs(cO).lt.eps)  goto  99 
al  =  0.5*(a0  +  bO) 
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bl  =  sqrt(aO*bO) 

cl  =0.5*(a0-b0) 

rmult  =  2.0*rmult 

stor  =  stor  +  rmult*(cl**2) 

bk=  1.5707963/al 

be  =  bk*(1.0  -  0.5*stor) 

aO  =  al 

bO  =  bl 

cO  =  cl 

goto  89 

♦ 

99  continue 
return 
end 

* 

Subroutine  Xsecenrg 

* 

*  This  subroutine  calculates  the  change  in  internal  energy  due  to  a 

*  change  in  the  cross-sectional  area  of  an  elliptic  tube  when  an 

*  internal  pressure  (greater  than  external  or  ambient  pressure)  is 

*  applied  to  the  tube.  The  change  in  internal  energy  occurs  when 

*  the  elliptical  cross-section  becomes  more  cicular  with  applied  pres- 

*  sure. 

* 

implicit  real*8  (a-h,o-z) 
real*8  mixs,  mib,  mu,  nu 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 

+  capl,capx,capy 

* 

call  jacconst 

♦ 

cape  =  (2.0*a*epr*mixs*cl*(c2dn  +  c3dn  +  c4dn))/(1.0  -  nu**2) 

* 

*  The  change  in  internal  energy  due  to  a  change  in  the  cross-section 

*  is  defined  as  duxs  =  capc*db**2 

* 

return 

end 
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Subroutine  Jacconst 


This  subroutine  simplifies  the  expression  needed  to  calculate  the 
change  in  internal  energy  due  to  a  change  in  the  tube's  cross- 
section.  Cl  is  a  constant  for  specific  major  and  minor  axes  of 
the  elliptical  tube.  Constants  c2nd,  c3dn,  and  c4dn  are 
also  specific  for  given  tube  axes  but  also  include  the  evaluation 
of  the  elliptic  trigonometric  function  dn(u)  raised  to  the  -4,  -6, 
and  -8  powers,  respectively. 

implicit  real*8  (a-h,o-z) 
real*8  mixs,  mib,  mu,  nu 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 

+  capl,capx,capy 

cl  =  (b/a)**4*(1.0/(be*a**2  -  bk*b**2)**2) 
c2  =  (((a**2  +  b**2)/b**2)*be  +  bk)**2 
c3  =  -6.0*be*(((a**2  +  b**2)/b**2)*be  +  bk) 
c4  =  9.0*be**2 

dn2  =  (2.0*be*(2.0  -  smk**2)/(1.0  -  smk**2)  -  bk)/(3.0*(1.0  - 
+  smk**2)) 

dn3  =  (4.0*(2.0  -  smk**2)/(15.0*(1.0  -  smk**2)**2))*(2.0*be*(2.0 
+  -  smk**2)/(1.0  -  smk**2)  -  bk)  -  (3.0*be/(5.0*(1.0  -  smk 

-I-  **2^**2)) 

dn4  =  (24.0*(2.0  -  smk**2)**2/(105.0*(1.0  -  smk**2)**3)  -  5.0/ 

+  (2 1 .0*(  1 .0  -  smk**2)**2))*(2.0*be*(2.0  -  smk**2)/(  1.0- 

+  smk**2)  -  bk)  -  ( 1 8.0*be*(2.0  -  smk**2)/(35.0*(  1 .0  -  smk**2)**3)) 

c2dn  =  c2*dn2 
c3dn  =  c3*dn3 
c4dn  =  c4*dn4 

return 

end 


Subroutine  Dummy 

This  subroutine  calculates  the  change  in  the  slope  at  the  free  end 
of  the  tube  using  Castigliano's  second  theorem.  The  change  in  the 
slope  is  determined  by  applying  a  dummy  moment  at  the  tube  tip  and 
incorporating  that  moment  into  the  internal  energy  equation  for  a 


beam.  Thetai  is  found  by  integrating  the  internal  energy  equation  and 
then  setting  the  dummy  moment  equal  to  zero. 

implicit  real  *8  (a-h,o-z) 
real  *8  mixs,  mib,nu,  mu 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda, da, dDdb,pi,dp,dR,db,dmu, thetai,  vofcapj , 

+  capl,capx,capy 

thetai  =  ((p  +  dp)*a*(R**2)*D*(mu  -  sin(mu)))/(epr*mib) 

Thetai  should  be  a  positive  value  to  indicate  that  the  angle  of 
rotation  has  the  same  sense  of  direction  as  the  applied  dummy  moment. 

return 

end 


Subroutine  Tip 

Once  thetai  is  known,  the  tip  deflection  of  the  tube  can  be 
calculated  by  assuming  that  mu*R  =  constant.  In  doing  so, 
mu*dR  +  R*dmu  =  0  must  be  true.  Geometry  shows  that  mu  has  the 
same  magnitude  but  opposite  sense  of  thetai. 

implicit  real  *8  (a-h,o-z) 
real  *8  mixs,  mib,  nu,  mu 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj , 

+  capl,capx,capy 

dmu  =  -thetai 

dR  =  (-R*dmu)/mu 

return 

end 


Subroutine  Bndenrg 


This  subroutine  simplies  the  expression  needed  to  calculate  the 
change  in  internal  energy  due  to  bending  in  the  tube.  Capa,  capb, 
and  capz  are  all  constant  for  the  given  parameters  but  change  as  the 
parameters  change. 

implicit  real*8  (a-h,o-z) 

real  *8  mixs,  mib,  nu,  mu 

common  a,b,smk,bk,be,epr,R,mixs,  mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 

+  capl,capx,capy 

capa  =  (((p  +  p)*(a**2)*(R**2)*D)*(6.0*musin(mu)+sin(2.0*mu)))/ 

+  (4.0*epr*mi) 


capb  =  (((p  +  dp)*R*b*(be  -  bk)/(be*a**2  -  bk*b**2))*(a*dDda  + 
+  D))  +  ((p  +  dp)*R*dDdb) 

capz  =  (R*D*dp)  +  (1.5*(p  +  dp)*D*dR) 

The  change  in  internal  energy  due  to  bending  is  defined  as 
dub  =  capa*(capb*db  +  capz) 

return 

end 


Subroutine  Dercapd 

This  subroutine  simplifies  the  expression  for  the  net  moment  acting 
on  the  tube  as  a  result  of  the  applied  internal  pressure.  D  is  the 
simplifying  term  in  the  moment  equation.  dDda  and  dDdb  are  the 
partial  derivatives  of  D  with  respect  to  a  and  b,  respectively. 

implicit  real*  8  (a-h,o-z) 
real  *8  mixs,  mib,  mu,  nu 

common  a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c  1 , 
+  c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 

+  capl,capx,capy 


D  =  a*((4.0/3.0)*((bk  -  be)/(smk**2)  +  (2.0*be  -  bk))  -  pi*(b/ 

+  a)) 

dDda  =  (4.0/3.0)*(((a**2*(bk  -  be)*(3.0*(a**2  -  b**2)  -  2.0*a** 
+  2))/(a**2  -  b**2)**2)  -  (2.0*be  -  bk)) 

dDdb  =  (4.0/3.0)*((2.0*a**3*b*(bk  -  be)/(a**2  -  b**2)**2)  -  pi) 

return 

end 
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